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ABSTRACT 

We  present  the  first  Spitzer  Infrared  Spectrograph  observations  of  the  [O  iv]  25.89  /xm  emission  line  detected  from 
the  ultraluminous  X-ray  source  (ULX)  in  Holmberg  IT  This  line  is  a  well-established  signature  of  high  excitation, 
usually  associated  with  active  galactic  nucleus  (AGN).  Its  detection  suggests  that  the  ULX  has  a  strong  impact 
on  the  surrounding  gas.  A  Spitzer  high-resolution  spectral  map  shows  that  the  [O  iv]  is  coincident  with  the  X-ray 
position  of  the  ULX.  The  ratios  of  the  [Oiv]  to  lower-ionization  lines  are  similar  to  those  observed  in  AGN, 
suggesting  that  a  strong  UV  and  X-ray  source  is  responsible  for  the  photoionization.  The  best  XMM-Newton  data 
are  used  to  model  the  X-ray  band  which  is  then  extrapolated  into  the  UV.  We  perform  infrared  and  ultraviolet 
photometry,  and  use  previously  published  optical  and  radio  data  to  construct  the  full  spectral  energy  distribution 
(SED)  for  the  ULX  and  its  companion.  The  preferred  model  to  describe  the  SED  includes  an  accretion  disk  which 
dominates  the  soft  X-rays  but  contributes  little  at  UV  and  optical  wavelengths.  The  optical  counterpart  is  consistent 
with  a  B  supergiant  as  previously  suggested  in  other  studies.  The  bolometric  luminosity  of  the  ULX  suggests 
the  presence  of  an  intermediate-mass  black  hole  with  mass  >85  Mq  for  sub-Eddington  accretion  or,  alternatively, 
a  stellar-mass  black  hole  that  is  accreting  at  super-Eddington  rates.  In  a  follow-up  second  paper,  we  perform 
detailed  photoionization  modeling  of  the  infrared  lines  in  order  to  constrain  the  bolometric  luminosity  of  the  ULX. 
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1.  INTRODUCTION 

The  majority  of  our  knowledge  of  ultraluminous  X-ray 
sources  (ULXs)  stems  from  over  200  observations  of  these 
objects  by  both  Chandra  and  XMM-Newton  (Liu  &  Mirabel 
2005).  The  elevated  X-ray  luminosities  of  ULXs  have  raised 
intriguing  questions  about  their  nature  and  the  source  of  the 
hard  X-ray  emission.  A  natural  explanation  for  the  substantial 
X-ray  luminosities  is  the  presence  of  intermediate-mass  black 
holes  (IMBH;  e.g.,  Colbert  &  Mushotzky  1999);  however, 
substantiating  evidence  for  such  objects  is  lacking  (Stobbart 
et  al.  2006;  Gon9alves  &  Soria  2006;  Berghea  et  al.  2008). 
Attention  has  switched  to  less  exotic  sources  for  the  emission, 
such  as  beaming  (King  et  al.  2001)  or  super-Eddington  accretion 
from  a  stellar-mass  black  hole  (Begelman  2002),  both  of  which 
explain  ULX  X-ray  properties  without  an  IMBH.  The  idea  is 
not  far  fetched,  since  examples  of  super-Eddington  sources  have 
been  found  in  our  Galaxy,  such  as  GRS  1915-^105  (Eender  & 
Belloni  2004),  V4641  Sgr  (Revnivtsev  et  al.  2002),  and  possibly 
SS  433.  Indeed  SS  433  could  be  an  example  of  both  beaming 
and  super-Eddington  accretion,  the  combination  of  which  could 
easily  explain  even  the  most  luminous  ULXs  (Begelman  et  al. 
2006;  Poutanen  et  al.  2007). 

Very  little  work  has  been  done  on  ULXs  at  non-X-ray 
wavelengths,  especially  in  comparison  with  active  galactic 
nucleus  (AGN)  and  Galactic  stellar-mass  black  holes.  This  is 
in  part  because,  unlike  standard  AGN,  the  emission  from  the 
accretion  process  itself  is  harder  to  detect  at  these  wavelengths. 
Indeed,  the  X-ray  to  optical  luminosity  ratios  in  ULXs  are  at 
least  10  times  higher  compared  to  those  seen  in  AGN  (Winter 
et  al.  2007).  To  make  matters  worse,  there  are  no  ULXs  in  our 
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own  Galaxy  (though  SS  433  is  regarded  as  a  ULX  by  some 
authors),  and  the  stellar  ULX  companions  are  very  difficult  to 
identify  in  other  galaxies. 

What  little  multi-wavelength  work  done  on  ULXs  thus 
far  has  been  executed  primarily  with  high-resolution  optical 
spectroscopy.  In  a  few  cases,  ULXs  have  been  found  to  have 
an  observable  effect  on  their  local  environment — developing 
ionized  bubble  nebulae  of  up  to  500  pc  diameter  (Pakull  & 
Mirioni  2002).  These  are  much  larger  than  regular  supernova 
remnants  (SNRs)  and  energy  estimates  show  that  the  bubbles 
are  consistent  with  being  blown  either  by  relativistic  jets  from 
a  ULX,  or  abnormally  powerful  supemovae.  In  some  cases, 
high-excitation  lines  have  been  detected,  such  as  Hen  A.4686 
(I.P.  =  54.4  eV;  e.g.,  Abolmasov  et  al.  2007).  Only  very  hot 
Wolf-Ray  et  stars  can  produce  this  line  (Schaerer  1996),  but 
the  ULX  companions  are  not  such  stars  (Abolmasov  et  al. 
2007).  Therefore,  the  detection  of  the  Hen  line  must  be  due 
to  the  interaction  of  the  ULX  with  the  surrounding  medium. 
Interestingly,  SS  433  is  also  surrounded  by  a  100  pc  radio  nebula, 
W  50  (e.g.,  Eabrika  2004).  Vazquez  et  al.  (2007)  conducted 
the  first  infrared  spectroscopy  study  of  ULXs.  They  obtained 
high-resolution  Spitzer  spectra  of  six  ULXs  in  NGC  4485/90. 
Vazquez  et  al.  (2007)  looked  for  the  mid-infrared  [Oiv]  high- 
ionization  line,  which  has  an  equivalent  ionization  potential 
to  He  II.  However,  apart  from  a  possible  weak  detection  in  one 
object  in  their  sample  (ULX-1 ,  which  is  probably  a  SNR),  [O  iv] 
was  not  found  in  the  NGC  4485/90  ULXs. 

The  importance  of  ionization  nebulae  for  investigating  X-ray 
sources  has  long  been  realized,  beginning  with  the  first  detection 
of  nebulae  around  supersoft  sources  such  as  LMC  X-1  (Pakull 
&  Angebault  1986)  and  CAL  83  (Remillard  et  al.  1995)  in  the 
Large  Magellanic  Cloud  (LMC).  The  high-excitation  lines  ob¬ 
served  in  nebulae,  such  as  the  He  ii  recombination  line,  can  act  as 
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giant  bolometers  for  the  ionizing  UV  and  X-ray  emission,  thus 
allowing  the  intrinsic  hard  spectral  energy  distribution  (SED) 
to  be  inferred.  Long-term  averaged  luminosities  can  be  esti¬ 
mated  and  the  UV  emission  constrained.  This  provides  tremen¬ 
dous  value  when  studying  the  nature  of  ambiguous  ULXs,  be¬ 
cause  the  luminosities  estimated  from  ionization  are  indepen¬ 
dent  of  X-ray  spectral  modeling  estimates.  In  fact,  such  lumi¬ 
nosities  are  more  accurate,  because  the  X-ray-based  estimates 
are  in  general  lower  bounds  and  depend  on  the  viewing  an¬ 
gle.  More  importantly,  the  ionization-estimated  luminosities  are 
time-averaged  estimates,  which  is  critical  for  variable  sources 
such  as  ULXs.  This  is  especially  true  for  ULXs  with  very  few 
X-ray  observations,  since  some  of  these  sources  can  be  tran¬ 
sient.  On  the  other  hand,  photoionization  modeling  is  depen¬ 
dent  on  many  parameters  which  are  often  not  known  very 
precisely,  such  as  the  gas  density,  geometry,  metallicity,  etc. 
Contamination  from  other  ionizing  sources  is  also  a  concern 
especially  in  crowded  regions.  However,  when  used  together. 
X-ray  analysis  and  photoionization  modeling  can  provide  good 
constraints  on  the  ULX  luminosity.  For  example,  a  simple  com¬ 
parison  between  the  two  estimates  can  constrain  the  beaming 
factor. 

One  of  the  most  interesting  examples  of  an  ionized  nebula  as¬ 
sociated  with  a  ULX  is  located  in  the  dwarf  galaxy  Holmberg  II 
(Pakull  &  Mirioni  2002;  Lehmann  et  al.  2005;  Abolmasov  et  al. 
2007).  Hen  is  detected  in  a  small  nebular  region  around  the 
ULX,  the  “Heel”  of  the  “Foot  Nebula”  (Kaaret  et  al.  2004, 
KWZ  hereafter).  An  OB  star  is  the  likely  counterpart  (KWZ). 
The  morphology  of  Hen  emission  and  other  optical  lines  is 
consistent  with  the  photoionization  of  the  nebula  by  the  ULX. 
KWZ  confirmed  this  with  CLOUDY  simulations.  Miller  et  al. 
(2005)  detected  extended  radio  emission  from  a  ~50  pc  diame¬ 
ter  region  coincident  with  the  position  of  the  ULX.  The  authors 
conclude  that  the  fiux  and  morphology  of  the  radio  emission  are 
inconsistent  with  emission  from  a  SNR  or  an  H  ii  region,  and  is 
probably  associated  with  the  ULX  activity. 

In  this  paper  and  a  subsequent  companion  paper,  we  use  IR 
Spitzer  spectroscopic  observations  for  the  first  time  to  constrain 
the  underlying  hard  SED  of  a  ULX.  The  IR  observations,  when 
combined  with  previously  published  optical  and  X-ray  data, 
suggest  that  the  powerful  UV  and  soft  X-ray  emission  from 
the  Holmberg  II  ULX  photoionizes  the  surrounding  medium.  In 
Section  2,  we  present  the  Spitzer  Infrared  Spectrograph  (IRS) 
data  analysis  andXMM-Newton  X-ray  spectral  fitting  results.  We 
perform  additional  photometry  measurements  with  the  Infrared 
Array  Camera  (IRAC),  the  Multiband  Imaging  Photometer  for 
Spitzer  (MIPS),  the  European  Photon  Imaging  Camera  (EPIC) 
and  Optical  Monitor  (OM)  instruments  on  XMM-Newton,  and 
the  Galaxy  Evolution  Explorer  {GALEX).  In  Section  3,  we 
analyze  the  possible  contribution  of  self-irradiated  accretion 
disk  to  the  UV  and  optical  fiuxes.  In  Section  4,  we  use  our  data 
and  other  published  data  to  construct  an  SED  of  the  Holmberg  II 
ULX.  This  is  also  used  in  the  companion  paper  to  further 
constrain  the  X-ray  spectrum.  In  Sections  5  and  6,  we  discuss 
the  significance  of  the  [O  iv]  line  detection,  compare  the  IR  line 
ratios  with  those  of  AGN,  and  present  the  final  observational 
results.  A  follow-up  paper  will  present  detailed  photoionization 
analysis  of  the  Spitzer  data. 

2.  OBSERVATIONS  AND  DATA  ANALYSIS 

We  begin  this  section  with  a  description  of  the  Spitzer 
IRS  (Houck  et  al.  2004)  observation  of  the  Holmberg  II 
ULX.  We  next  present  multi-wavelength  photometry  and 


X-ray  spectroscopy  used  in  the  subsequent  section  to  construct 
the  SED  for  the  ULX. 

2.7.  Spitzer  IRS  Observations 

The  data  presented  here  are  archival  Spitzer  data,  taken  on 
2004  October  4,  as  part  of  the  Legacy  Program  Spitzer  Infrared 
Nearby  Galaxies  Survey  (SINGS  Program  193;  see  Kennicutt 
et  al.  2003).  The  IRS  spectral  map  is  described  by  three  parallel 
pointings  and  five  perpendicular  at  half-slit  steps,  for  both 
Short  High  (SH,  in  the  range  9.9-19.6  /xm)  and  Long  High 
(LH,  in  the  range  18.7-37.2  /xm)  modules.  The  total  exposure 
time  was  240  s  for  both  SH  and  LH  observations.  The  data  were 
preprocessed  by  the  IRS  pipeline  (ver.  1 1)  at  the  Spitzer  Science 
Center  prior  to  download.  Spectral  maps  were  constructed 
using  “BCD-level”  processed  data  in  conjunction  with  CUBISM 
(cube  Builder  for  IRS  Spectra  Maps)  version  1.5^  (Smith  et  al. 
2007;  Kennicutt  et  al.  2003).  The  absolute  calibration  error  on 
fiuxes  obtained  with  CUBISM  is  ~25%  (Smith  et  al.  2004). 

We  adopt  a  distance  of  3.05  Mpc  to  Holmberg  II  (KWZ;  V' 
is  15  pc).  The  LH  maps  obtained  with  CUBISM  have  intrinsic 
sizes  of  44('6  x  28('8  (670  pc  x  430  pc),  while  the  SH  maps 
are  21"  x  13('6  (405  pc  x  205  pc).  However,  in  order  to  avoid 
aperture  effects,  the  spectra  presented  here  were  extracted  from 
like-sized  apertures.  Thus,  the  SH  and  LH  spectral  data  result 
from  an  equally  sized  aperture  of  23''  x  13'.'6  (340  pc  x  205  pc), 
which  is  only  slightly  smaller  than  the  full  SH  aperture.  We 
also  extract  spectra  from  the  smallest  aperture  permitted  by 
CUBISM:  an  8'.'9  x  8'.'9  aperture  centered  on  the  ULX  (i.e.,  2  x 
2  LH  pixels,  or  1 30  pc  x  1 30  pc).  We  will  refer  to  these  apertures 
as  the  standard  aperture  and  the  small  aperture,  respectively. 

The  high-resolution  spectra  (10-37  /xm)  for  the  two  apertures 
are  plotted  in  Figure  1.  Since  the  spectra  are  extracted  from 
matched  apertures,  the  spectral  lines  in  both  the  SH  and  LH 
modules  are  properly  calibrated  and  the  continuum  is  properly 
matched.  In  Figure  2,  we  show  in  detail  the  four  detected 
lines,  with  the  theoretical  line  centers.  We  used  SMART 
(Spectroscopic  Modeling  Analysis  and  Reduction  Tool;  see 
Higdon  et  al.  2004)  version  6.2.6  to  measure  line  fiuxes.  Only 
four  lines  were  detected,  the  fiuxes  for  these  and  upper  limits 
for  nondetections  are  presented  in  Table  1 .  Here,  we  show  the 
results  for  both  the  standard  aperture  and  the  small  aperture. 
We  note  that  the  [O  iv]  line  is  detected  at  a  signal-to-noise  ratio 
(S/N)  >10.  The  spectral  lines  for  the  Holmberg  II  galaxy  are 
also  presented  in  Dale  et  al.  (2009).  The  focus  of  that  paper 
was  not  the  ULX  in  Holmberg  II,  and  as  a  result,  the  extraction 
region  for  the  relevant  spectral  lines  is  larger  than  the  apertures 
presented  here.  However,  with  this  in  mind,  we  find  that  our  line 
fiuxes  agree  reasonably  well  with  those  presented  in  Dale  et  al. 
(2009). 

The  SH  and  LH  maps  for  the  four  detected  lines  are  shown 
in  Figure  3.  A  false-color  Hubble  Space  Telescope  (HST)  image 
shows  the  Foot  Nebula  in  Ha  emission,  together  with  the 
standard  and  small  extraction  apertures.  Because  of  the  poor 
spatial  resolution  of  the  Spitzer  spectral  maps,  the  emission 
for  most  lines  is  difficult  to  disentangle  from  the  nearby  star¬ 
forming  regions  in  the  Foot  Nebula.  This  is  especially  true 
of  the  low-ionization  lines  ([Neiii],  [Siii],  and  [Sin]).  The 
morphologies  of  these  maps  indicate  that  the  emission  is  likely 
coming  from  both  the  ULX  and  the  surrounding  star-forming 
regions.  However,  while  the  low-ionization  lines  [Sin]  and 
[S  III]  are  very  extended,  the  [O  iv]  emission  is  visibly  compact. 


^  Available  at:  http://ssc.spitzer.caltech.edu/archanaly/contributed/cubism/ 
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Figure  1.  Full  IRS  spectra  obtained  from  matched  apertures  with  CUBISM  from  the  two  apertures  shown  in  Figure  3.  (a)  The  standard  aperture,  just  slightly  shorter 
than  the  SH  aperture,  (b)  The  small  aperture,  this  is  the  smallest  extraction  aperture  possible  for  matched  spectra,  4  LH  pixels.  By  comparing  the  spectra,  it  is  evident 
that  the  [O  iv]  is  clearly  concentrated  in  the  small  aperture. 


Wavelength  (^im)  Wavelength  (ftm) 


Figure  2.  Details  of  the  detected  infrared  lines  in  the  spectrum  of  Figure  1(a).  The  vertical  lines  show  the  theoretical  center  for  the  lines  at  rest  (dotted)  and  redshifted 
(solid).  We  used  the  recessional  velocity  adopted  by  KWZ  (157  km  s“^). 


The  peak  of  the  of  the  [Oiv]  emission  is  coincident  with  the 
Chandra  position  for  the  ULX,  with  the  majority  of  the  emission 
concentrated  in  the  2  x  2  pixel  region,  overlapping  with  the 
X-ray  point  source.  There  is  a  slight  extension  to  the  west  and 
south  that  is  visible  on  the  spectral  map  and  which  may  imply 
that  the  [Oiv]  emission  is  extended.  These  faint  extensions 
correspond  very  well  with  two  regions  of  young  star  formation 
visible  in  the  HST  image  in  Eigure  3. 

KWZ  found  that  the  Hen  line,  which  has  an  ionization 
potential  just  slightly  lower  than  [Oiv]  (54.4  eV,  compared 
to  54.9  eV),  is  also  concentrated  in  a  small  region  of  50  pc 
diameter  around  the  identified  optical  counterpart  of  the  ULX. 
Because  the  spatial  resolution  is  limited,  a  detailed  comparison 
of  our  maps  with  the  optical  maps  presented  by  KWZ  cannot  be 
made.  However,  we  find  that  the  morphology  of  our  [O  iv]  map 
appears  similar  to  He  ii,  being  brighter  in  the  east  and  fainter 
and  more  extended  in  the  west. 


2.2.  X-ray  Observations  and  Spectral  Fitting 

We  use  X-ray  spectroscopy  to  model  the  X-ray  emission 
from  the  ULX.  We  also  use  it  to  estimate  the  UV  emission  by 
extrapolating  the  model  to  lower  energies.  X-ray  observations 
of  the  Holmberg  II  galaxy  are  available  from  ROSAT,  ASCA, 
XMM-Newton,  and  Chandra  (Zezas  et  al.  1999;  Miyaji  et  al. 
2001;  Kerp  et  al.  2002;  Kaaret  et  al.  2004;  Goad  et  al.  2006). 
The  X-ray  luminosity  of  the  ULX  is  in  excess  of  10"^^  erg  s“^  in 
all  observations,  with  the  exception  of  a  peculiar  low  and  soft 
state  in  2002  September,  when  the  luminosity  dropped  to  ^^5  x 
1039  gj.g  g-i  (Dewangan  et  al.  2004). 

To  model  the  X-ray  spectrum,  we  used  the  best  available 
data  from  the  longest  (~100  ks)  XMM-Newton  EPIC  PN 
observation  in  the  archive  (Obs.  ID  0200470101,  taken  on 
2003  December  5).  The  data  were  processed  by  Winter  et  al. 
(2006).  We  used  XSPEC  version  12.4  to  model  the  spectrum 
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Figure  3.  HST  ACS  false-color  image  of  the  Foot  Nebula  in  Holmberg  II  (left),  and  the  Spitzer  spectral  maps  of  the  infrared  detected  lines  (right).  The  colors  in  the 
HST  image  correspond  to  the  following  hlters:  814W  (red),  658N  (green.  Ho;),  and  550M  (blue).  The  spectral  maps  are  [Neiii]  15.56  /xm,  [Siii]  18.71  /xm,  [Oiv] 
25.89  /xm,  [S  iii]  33.48  /xm,  and  [Si  ii]  34.82  /xm.  The  units  on  the  color  bar  are  W  m~^  sr“^  The  apertures  used  to  extract  IRS  spectra  (standard  and  the  small  aperture, 
see  Section  2.1)  are  overplotted  on  all  images.  In  the  HST  image,  we  also  show  the  whole  LH  map  for  comparison  with  the  LH  maps  shown  at  right.  The  ULX  optical 
counterpart  identified  by  KWZ  has  the  appearance  of  a  bright  star  in  the  HST  image,  close  to  the  center  of  the  small  aperture,  and  its  position  is  shown  with  a  green 
X  on  the  spectral  maps.  It  is  located  at  the  center  of  the  “Heel” — the  round-shaped  strong  Ho;  emission.  An  OB  association  of  stars  is  seen  west  from  the  ULX  (see 
also  the  OM  image  in  Figure  4).  They  appear  to  have  blown  a  large  bubble  in  the  gas,  creating  the  “Sole”  and  “Toes”  of  the  Foot  Nebula.  The  low-ionization  lines  are 
probably  contaminated  by  these  stars.  Some  contamination  is  also  produced  by  the  bright  foreground  star,  seen  on  the  HST  image  just  inside  the  standard  aperture, 
close  to  its  northern  edge. 


Table  1 

Measured  Infrared  Lines 


Line 

Standard  Aperture 

Small  Aperture 

Flux 

S/N 

L 

Flux 

S/N 

L 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

[Neii]  12.81  /xm 

<6.08 

<0.68 

<1.98 

<0.22 

[Neiii]  15.56  /xm 

12.15  ±3.4 

4.60 

1.35  ±0.39 

5.46  ±  1.37 

7.6 

0.61  ±0.15 

[Sill]  18.71  /xm 

2.64  ±  1.14 

4.64 

0.29  ±0.13 

3.02  ±  1.16 

4.7 

0.34  ±0.13 

[Nev]  24.32  /xm 

<3.88 

<0.43 

<1.16 

<0.13 

[Oiv]  25.89  /xm 

7.01  ±  1.60 

10.16 

0.78  ±  0.2 

3.63  ±  0.91 

44.4 

0.40  ±0.1 

[Fell]  25.99  /xm 

<2.13 

<0.24 

<0.78 

<0.09 

[Sill]  33.48  /xm 

<4.39 

<0.49 

<1.53 

<0.17 

[Sin]  34.82  /xm 

21.91  ±5.48 

9.90 

2.44  ±  0.61 

5.87  ±  1.47 

7.2 

0.65  ±0.16 

Notes.  We  show  line  fluxes  for  the  standard  (slightly  smaller  than  the  SH  map)  and  the  small  aperture  (4  LH  pixels)  as  defined  in 
Section  2.  Fluxes  are  in  10“^^  W  cm“^,  luminosities  (L)  in  10^^  erg  s“^ .  If  the  measurements  errors  are  smaller  than  the  absolute 
calibration  accuracy  of  25%,  the  latter  were  used.  For  nondetections  we  show  3a  upper  limits. 


between  0.3  and  8.0  keV.  The  X-ray  spectrum  is  the  most 
important  component  of  the  total  SED  since  it  is  the  only 
observed  emission  that  directly  emanates  from  the  ULX,  and 
is  not  contaminated  by  the  surrounding  star-forming  region.  For 
this  reason,  we  have  taken  many  precautions  in  modeling  the 
X-ray  spectrum. 

The  (far)  UV  emission  from  the  ULX  is  the  least  known,  but 
the  most  important  for  ionization.  The  only  available  method 
to  estimate  the  shape  of  the  UV  is  to  extrapolate  from  X-ray 
models.  Therefore,  it  is  very  important  to  find  the  right  X-ray 
model  to  fit  the  data,  especially  in  the  soft  X-ray  band.  Three 
physically  different  models  were  used.  Each  model  contains 
different  information  about  the  accretion  process.  The  first 
model  is  very  basic  and  is  widely  used  to  model  the  accretion 
process  for  stellar-mass  black  holes,  ULXs,  and  other  accreting 


objects  (Miller  et  al.  2004;  McClintock  &  Remillard  2006; 
Winter  et  al.  2006;  Feng  &  Kaaret  2006).  The  second  is  a 
modification  of  the  first  model  that  is  more  physically  realistic 
when  extrapolated  to  lower  energies.  Finally,  the  third  model 
is  more  flexible  than  the  other  two  models  in  the  UV  and 
allows  for  a  tighter  fit  in  the  UV  while  staying  within  the 
X-ray  error  bars.  The  intrinsic  (unabsorbed)  X-ray  luminosity  is 
>  10"^^  erg  s“^  for  all  models.  The  results  from  the  X-ray  spectral 
fits  are  listed  in  Table  2,  and  a  summary  of  the  models  is  as 
follows: 

1.  Power-law  Plus  Multicolor  Disk  (PLMCD)  model.  This 
is  a  two-component  model  with  a  multi-color  accretion 
disk  component  (MCD;  Mitsuda,  et  al.  1984)  and  a  power 
law  for  the  hard  tail.  Following  KWZ,  the  absorption 
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Table  2 

X-ray  Model  Fits 


Model 

(1) 

A^h 

(2) 

(10^^  cm“^) 

kTin/Fi 

(3) 

(keV) 

r/r2 

(4) 

Ax^/dof 

(5) 

L 

(6) 

(10^^  erg  s“^) 

MCD  fiux 

(7) 

Modified  PLMCD 

3.26  ±  0.56 

0.38  ±  0.02 

2.54  ±  0.02 

1.02/135,0.974/713 

1.11 

0.49 

PLMCD 

10.4  ±  0.70 

0.27  ±  0.02 

2.42  ±  0.05 

1.1386/849 

2.28 

0.12 

BPL 

5.49  ±  0.78 

1.73  ±0.10 

2.57  ±  0.02 

1.0316/849 

1.41 

0.60 

Notes.  (1)  X-ray  model,  as  described  in  Section  2.2.  (2)  Intrinsic  hydrogen  column  density.  The  Galactic  column  density  from  KWZ  (3.42  x  10^^  cm“^) 
was  added  separately.  (3)  Inner  disk  temperature  for  the  MCD  component.  For  the  BPL  model,  this  is  the  first  (low  energies)  photon  index  parameter. 
The  break  is  at  1 .0  keV.  (4)  Photon  index  for  the  power-law  component.  For  the  BPL  model,  this  is  the  second  (high  energies)  photon  index  parameter. 
(5)  Reduced  y}  values  for  the  fit  and  the  number  of  degrees  of  freedom.  For  the  Modified  PLMCD  model,  we  quote  two  values  for  each  separately 
htted  component  (see  the  text  for  details).  (6)  Unabsorbed  (intrinsic)  luminosities  between  0.1  and  10  keV.  (7)  MCD  component  unabsorbed  flux  as 
fraction  of  the  total  flux  for  the  PLMCD  model.  For  the  other  two  models,  this  simply  the  ratio  of  the  flux  between  0.1  and  1.0  keV  to  the  flux  in  the 
whole  range  (0.1-10  keV). 


column  density  was  fixed  at  the  Galactic  value  of  3.42  x 
10^^  cm“^,  and  an  additional  local  absorption  (Ah)  model 
was  added.  We  assume  a  local  metallicity  of  0.1  Zq  (see 
Kaaret  et  al.  2004;  Pilyugin  et  al.  2004).  This  is  a  very 
simplistic  model  commonly  used  to  fit  X-ray  data  for 
ULXs.  However,  the  model  falls  apart  when  extrapolated 
to  lower  energies,  where  the  power  law  dominates.  The 
extrapolation  is  actually  physically  unrealistic,  since  it  is 
impossible  for  the  disk  component  to  Comptonize  lower 
energy  photons.  Done  &  Kubota  (2006)  pointed  out  that  the 
power-law  component  in  such  cases  needs  to  have  a  break 
at  low  energies,  in  agreement  with  the  theoretical  model 
where  the  photons  from  the  disk  are  upscattered  in  the  hot 
corona.  Moreover,  if  the  power  law  is  produced  by  another 
mechanism  (e.g.,  a  jet),  it  would  not  extend  down  to  the  soft 
UV  and  optical  spectrum  because  the  fluxes  here  would  be 
over  predicted  (see  below.  Section  4).  Nevertheless,  the 
model  continues  to  be  used  extensively  in  X-ray  spectral 
fitting  of  ULXs.  In  conclusion,  this  model  does  fit  the  X-ray 
data  well  (Ax^/dof  =  1.14),  but  it  does  not  make  physical 
sense  when  extrapolated  to  lower  energies.  It  is  shown  here 
as  a  base,  from  which  alterations  can  be  made,  resulting  in 
more  physically  meaningful  models.  The  following  model 
is  one  such  modification  to  the  PLMCD  model,  and  is  more 
physically  realistic  when  extrapolated  into  the  UV. 

2.  Modified  PLMCD  model.  In  order  to  isolate  the  intrinsic 
disk  emission  from  the  Comptonized  disk  emission,  we 
fitted  the  MCD  and  power  law  separately  for  two  different 
regions  of  the  X-ray  spectrum.  We  first  fitted  the  hard  tail 
with  a  power  law.  The  spectrum  departs  from  this  power  law 
around  ~1  keV.  Then,  we  fitted  an  MCD  model  to  the  0.3- 
1.0  keV  band.  The  disk  temperature  (0.38  keV)  is  slightly 
higher  than  other  published  fits  (KWZ  obtained  ~0.2  keV), 
because  our  fitting  methods  are  slightly  different.  Lower 
temperature  accretion  disks  are  generally  more  effective 
ionizing  sources.  Therefore,  this  higher  temperature  model 
may  be  viewed  as  a  lower  bound  to  the  SED  in  the  EUV  and 
soft-X-ray  bands.  We  note  that  the  Modified  PLMCD  differs 
from  the  PLMCD  model  in  that  the  UV  band  is  predicted 
from  the  disk  component  rather  than  the  extrapolated  PL 
component.  This  difference  yields  a  much  more  physically 
realistic  interpretation  for  the  emission. 

3.  Broken  power-law  (BPL)  model.  Linally,  we  also  fitted  the 
X-ray  spectrum  to  a  simple  BPL  model.  Kaaret  et  al.  (2003) 
fit  a  similar  model  to  the  ULX  in  NGC  5408.  Such  a  model 
could  be  explained  if  the  emission  from  the  companion  star 


were  Comptonized  by  a  jet  (Georganopoulos  et  al.  2002). 
This  model  is  much  more  flexible  in  the  UV  than  either 
the  PLMCD  or  the  Modified  PLMCD  models,  and  can 
predict  the  slope  of  the  SED  in  the  UV  while  simultaneously 
fitting  the  hard  power-law  component  in  the  X-rays  within 
the  error  bars  of  the  X-ray  measurements.  Surprisingly, 
this  simple  model  provides  a  very  good  fit  to  the  X-ray 
data,  better  than  the  PLMCD  model.  The  break  energy  is 
at  1.0  keV,  the  same  energy  at  which  we  split  the  two 
components  in  the  Modified  PLMCD  model  (2  above). 
Linally,  we  note  that,  even  if  this  model  has  a  power  law 
less  steep  than  the  PLMCD  model  at  lower  energies,  it 
is  still  overpredicting  the  soft  UV  and  optical  emission 
according  to  our  photometry  results  (Section  4).  A  similar 
overprediction  at  lower  energies  was  found  for  this  model 
by  Done  &  Kubota  (2006)  for  the  Galactic  binary  XTE 
J1550— 564  when  at  high  accretion  rates  (in  the  very  high 
state). 

2.3.  SED  of  the  ULX:  Multi-wavelength  Photometry 

In  order  to  constrain  ULX  multi-wavelength  SED,  we  per¬ 
formed  photometry  in  the  IR  and  UV  using  Spitzer,  XMM- 
Newton,  and  GALEX  data.  We  also  incorporated  previously  pub¬ 
lished  optical  and  radio  data  from  HST  and  the  Very  Large  Array 
(VLA).  A  short  description  of  the  data  and  the  measurements 
are  shown  in  Table  3.  All  our  photometric  measurements  were 
extracted  from  an  aperture  of  50  pc  diameter  (3'.'3).  We  used  the 
optical  counterpart  identified  by  KWZ  to  center  the  extraction 
aperture  (see  Ligure  3).  Aspect  correction  for  the  IR  and  UV 
images  was  performed  using  bright  sources  in  field,  such  as  the 
foreground  star  seen  in  Ligure  3.  The  measurements  presented  in 
Table  3  represent  the  integrated  fluxes  over  this  region,  including 
the  nebular  emission  and  the  contamination  from  other  nearby 
sources.  They  should  therefore  be  considered  upper  bounds  to 
the  intrinsic  ULX  contribution. 

IR  photometry.  IRAC  and  MIPS  images  (see  Ligure  4)  were 
taken  as  a  part  of  SINGS  Program  195.  We  extracted  3.6  /xm 
and  8.0  /xm  fiuxes  using  50  pc  diameter  apertures  centered  on 
the  ULX  position.  Lor  the  MIPS  24  /xm  photometry,  we  used 
a  100  pc  extraction  aperture  and  then  divided  the  flux  by  four 
to  obtain  an  approximate  fiux  for  an  equivalent  50  pc  aperture. 
We  followed  the  standard  extraction  procedure  and  conversion 
factors  provided  in  the  IRAC  Data  Handbook.^  The  pipeline 
mosaics  for  IRAC  and  MIPS  are  accurate  to  within  20%. 


^  http :// s sc .  spitzer. caltech . edu/irac/ dh/ 
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Table  3 

Multi- wavelength  Photometry  Measurements 


Instrument 

Program  (Observation) 

Wavelength 

Date 

Flux  Density 

Corrected  Flux 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Photometry  performed  in  this  paper 

Spitzer  IRAC 

SINGS  159 

3.6  /xm 

2004  Oct  10 

0.15  ±  0.03  mJy 

Spitzer  IRAC 

SINGS  159 

8.0  /xm 

2004  Oct  10 

0.26  ±  0.05  mJy 

Spitzer  MIPS 

SINGS  159 

24.0  /xm 

2004  Oct  14 

0.69  ±  0.14  mJy 

XMM-Newton  OM 

0200470101 

UVWl  2910  A 

2004  Apr  15 

33.1  ±  0.1  /xJy 

0.60 

57.5  /xJy 

NUV 

NGS 

1750-3013  A 

2003  Dec  5 

26.3  ±  0.2  /xJy 

0.95 

63.1  /xJy 

GALEX  FUV 

NGS 

1350-2013  A 

2003  Dec  5 

33.7  ±  0.3  /xJy 

0.84 

73.0  /xJy 

Previously  published  photometry 

VLA 

AT159 

1.4  GHz 

1994  Apr  3 

1.174  ±  0.085  mJy 

VLA 

AR165,  AR227,  AR258 

4.86  GHz 

1990  Jul  21-1991  Dec  6 

0.667  ±  0.207  mJy 

HST  ACS 

GO  9684 

F550M 

2002  Nov  24 

6.56  ±  0.03  /xJy 

0.32 

9.51  /xJy 

Notes.  All  the  photometric  measurements  performed  in  this  paper  used  an  aperture  of  50  pc  diameter  (3"3).  These  values  represent  the  integrated 
fluxes  over  this  region,  including  the  nebular  emission  and  are  contaminated  by  other  nearby  sources,  therefore  we  consider  them  as  upper  bounds.  The 
published  data  are  taken  from  Miller  et  al.  (2005)  (radio  VLA)  and  from  KVZ  (optical  HST).  The  radio  emission  is  extended,  and  is  comparable  to  our 
extraction  region.  The  optical  measurement  is  for  the  star-like  source  which  is  believed  to  be  the  companion  of  the  accreting  BH.  This  measurement  is 
likely  less  contaminated  by  the  nebular  emission  or  other  sources.  The  extinctions  in  Column  6  are  both  Galactic  and  local  from  KWZ.  For  the  infrared 
an  radio  data,  the  extinctions  are  negligible. 


UV  photometry.  XMM-Newton  OM  images  of  Holmberg  II 
(Figure  4)  are  available  in  the  UVW 1  band,  which  is  centered  at 
2910  A  (FWHM  620  A).  We  use  the  mosaicked  image  created 
for  the  Optical  Monitor  Catalog  (OMCat;  Kuntz  et  al.  2008), 
which  has  an  improved  coordinate  correction  (better  than  O'.' 5) 
over  the  automated  pipeline.  To  measure  the  UVW  1  count  rate, 
we  used  the  XIMAGE  software  included  in  the  HE  As  oft  package. 
We  converted  the  count  rate  into  flux  using  the  conversion 
factors  listed  in  the  XMM-Newton  Science  Analysis  System 
User’s  Guide f  We  obtained  a  UVWl  magnitude  of  18.42  zb 
0.06,  which  corresponds  to  a  luminosity  of  8.3  X  10^^  erg  s  ^ 
and  a  flux  density  of  33.1  zb  0.1  /xJy.  The  OM  image  in  Figure  4 
shows  that  the  Foot  Nebula  has  a  strong  UV  emission  which 
peaks  on  the  OB  association  of  stars  seen  in  the  HST  image  of 
Figure  3  west  of  the  UFX.  The  total  UVWl  flux  for  this  UV 
bright  region  (including  the  UFX  and  its  optical  counterpart)  is 
equivalent  to  29  05V  stars.  Figure  3  shows  that  these  bright  stars 
have  apparently  blown  a  cavity  in  the  NW  part  of  the  nebula, 
thus  producing  the  “Sole”  and  the  “Toes”  of  the  Foot  Nebula. 
We  note  that  these  cavities  are  correlated  with  the  faint  extended 
[Oiv]  emission  seen  in  the  CUBISM  spectral  maps  (Figure  3). 
The  strongest  emission  of  Ha  is,  however,  in  the  Heel  that  hosts 
the  UFX. 

We  also  used  archived  GALEX  observations,  taken  as  a  part 
of  the  GALEX  Nearby  Galaxies  Survey  (NGS;  Bianchi  et  al. 
2003).  Photometric  fluxes  were  extracted  from  the  FUV  (1350- 
20 1 3  A)  and  NUV  ( 1 750-20 1 3  A)  bands  with  the  XIMAGE  tool. 
The  GALEX  images  have  a  lower  resolution  compared  to  OM, 
but  they  show  similar  features. 

Published  radio  data.  Finally,  we  searched  the  literature 
for  published  photometric  data  for  the  Holmberg  II  UFX. 
Using  radio  VFA  data.  Miller  et  al.  (2005)  measured  1.174  zb 
0.0085  mJy  at  1.4  GHz  and  0.677  zb  0.207  mJy  at  4.86  GHz. 
Miller  et  al.  (2005)  found  extended  emission  covering  the 
Chandra  position  of  the  UFX:  3'.'7  x  2'.'7  at  1.4  GHz  and  3'.'4  x 
1'.'9  at  4.86  GHz.  These  are  comparable  to  the  50  pc  extraction 
region  used  in  the  IR  and  UV  photometry. 


http  ://xmm2 .  esac .  esa.int/ external/xmm_user_support/ documentation/ 


Figure  4.  Spitzer  IRAC  (8  /xm  Alter,  left)  and  XMM-Newton  OM  (UVW  1  Alter, 
right)  images  of  the  Foot  Nebula.  The  aperture  used  for  the  ULX  photometry 
is  shown  as  a  small  circle.  At  the  Holmberg  II  distance,  it  has  a  diameter  of 
50  pc.  With  dotted  lines,  we  show  the  two  apertures  used  to  extract  the  IRS 
spectra  shown  in  Figure  I .  The  larger  aperture  (continuous  line)  corresponds  to 
the  whole  LH  map.  Notice  that  the  bright  star  which  contaminates  the  infrared 
emission  in  Figure  3  is  very  faint  in  the  UV  image.  The  strongest  UV  emission  in 
the  nebula  is  located  west  of  the  ULX  and  it  originates  from  the  OB  association 
better  seen  in  Figure  3. 


Published  optical  photometry.  In  addition,  KWZ  measured 
a  V-band  magnitude  of  22.04  zb  0.08  for  the  stellar  optical 
counterpart  from  observations  with  the  Advanced  Camera  for 
Surveys  (ACS;  F550M  Alter,  see  Figure  3)  on  the  HST.  The 
optical  flux  is  likely  less  contaminated  by  the  nebular  emission 
or  emission  from  other  sources  compared  to  our  measurements, 
because  of  the  much  smaller  aperture. 

Extinction  corrections.  The  Milky  Way  reddening  in  the 
direction  of  Holmberg  II  is  E{B  —  V)  =  0.032  (from  the  NASA 
Extragalactic  Database  (NED)).  KWZ  estimated  an  additional 
local  reddening  of  E(B  —  V)  =  0.07  zb  0.01  based  on  the 
X-ray  absorption.  Thus,  the  total  extinction  for  Ry  =  3.1 
is  Ay  =  0.32.  We  show  the  projected  extinction  and  the 
extinction-corrected  values  for  the  UV  and  optical  data  in 
the  last  two  columns  of  Table  3.  These  were  calculated  using 
the  Cardelli  et  al.  (1989)  extinction  curves.  For  the  infrared  data, 
we  estimate  from  Draine  (1989)  that  the  extinction  is  Ax  <  0.02, 
and  therefore  negligible.  We  note  that  the  photometric  fluxes 
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Figure  5.  Examples  of  irradiated  disks  spectra.  The  higher  flux  spectrum  is 
calculated  for  a  IOOOMq  intermediate-mass  BH  binary  and  an  outer  radius 
of  10^^  cm.  The  lower  spectrum  is  for  a  stellar- mass  BH  binary  of  IOMq 
and  two  disk  sizes:  3  x  10^^  cm  and  10^^  cm.  The  continuous  lines  show 
the  total  spectrum,  dotted  lines  the  spectrum  from  inner  disk,  not  affected 
by  the  irradiation,  and  the  dashed  lines  from  the  outer  part.  The  vertical  line 
marks  the  [O  iv]  ionization  edge.  The  reprocessed  photons  have  energies  below 
10  eV,  and  therefore  not  ionizing.  All  cases  assume  accretion  at  Eddington 
luminosity. 


presented  here  are  likely  contaminated  by  the  surrounding 
nebular  emission  as  is  evidenced  by  the  HST  images  in  Eigure  3. 
Moreover,  KWZ  caution  that  the  local  extinction  could  in  fact 
occur  very  close  to  the  compact  object  and  therefore  not  affect 
our  measurements  at  all.  We  will  therefore  take  the  (extinction- 
corrected)  photometric  measurements  as  upper  bounds. 

3.  SELE-IRRADIATION  OE  THE  ACCRETION  DISK 

While  the  MCD  accretion  disk  model  is  a  good  approximation 
of  the  UV  continuum  to  first  order,  it  does  not  account  for 
self-irradiation  of  X-ray  photons  caused  by  thickening  of  the 
accretion  disk  at  the  outer  edges  and  relativistic  effects  on  the 
X-rays  emitted  from  the  inner  regions  of  the  disk.  In  Galactic 
X-ray  binaries,  this  refiected  radiation  has  been  shown  to 
increase  the  optical  luminosity  of  the  accretion  disk  by  as  much 
as  18%  (Sanbuichi  et  al.  1993;  Vrtilek  et  al.  1990;  Revnivtsev 
et  al.  2002).  Thus,  if  irradiation  were  contaminating  the  UV 
emission  in  the  Holmberg  II  ULX,  such  a  model  must  be 
included  in  the  final  fit. 

We  therefore  checked  to  see  if  X-ray  irradiation  would  affect 
the  ionizing  UV  emission  from  the  Holmberg  II  ULX,  or 
contribute  significantly  to  the  fiux  of  the  optical  counterpart. 
Eollowing  Sanbuichi  et  al.  (1993),  we  modeled  the  expected 
SED  of  an  irradiated  disk  for  black  holes  of  mass  10,  100, 
and  IOOOMq  and  with  varying  accretion  disk  radii.  Here,  the 
accretion  disk  radius  is  constrained  by  the  size  of  the  Roche  lobe, 
and  is  therefore  dependent  on  the  spectral  type  of  the  companion 
star.  We  find  that  for  stellar-mass  black  holes,  irradiation  affects 
the  accretion  disk  model  significantly  for  most  companion 
spectral  types.  However,  for  intermediate  black  hole  masses, 
irradiation  is  only  important  to  systems  with  late-type  supergiant 
companions. 


In  Eigure  5,  two  examples  of  irradiated  spectra  are  plotted. 
These  represent  typical  spectra  for  stellar-mass  BHs  (IOMq) 
and  IMBHs  (IOOOMq).  In  the  latter  case,  an  accretion  disk 
typical  for  a  system  with  an  M-supergiant  was  used  to  illustrate 
the  effect  of  irradiation.  Only  if  the  IMBH  has  such  companions 
is  the  self-irradiation  effect  significant.  However,  we  note  that  an 
M-supergiant  companion  is  not  consistent  with  the  observations 
for  the  ULX  in  Holmberg  II  (see  Section  4).  As  can  be  seen  from 
this  plot,  in  both  cases  the  irradiation  dominates  the  optical 
and  near-IR  spectrum  much  more  than  the  UV  band.  Thus, 
we  conclude  that  while  irradiation  may  be  important  to  model 
in  the  optical  spectrum  of  a  system  with  a  stellar-mass  black 
hole,  it  is  likely  not  dominant  for  an  IMBH  system.  Einally, 
we  conclude  that  the  irradiation  does  not  affect  the  detected  IR 
lines.  Moreover,  because  the  measured  fiuxes  are  much  higher 
than  the  predictions  from  the  accretion  disk,  self-irradiation  is 
not  likely  to  make  a  significant  contribution  to  our  measured 
fiuxes  at  any  wavelength. 

4.  SPECTRAL  ENERGY  DISTRIBUTION  (SED) 

Using  the  measurements  and  the  X-ray  models  in  Section  2, 
we  construct  the  SED  of  the  ULX  using  the  IR,  optical, 
UV,  and  X-ray  photometric  measurements.  The  final  SED 
is  plotted  with  the  continuous  black  line  in  Eigure  6.  All 
photometric  measurements  described  in  the  previous  section 
were  extracted  from  an  aperture  of  ~50  pc  centered  on  the 
ULX  (Eigure  4).  The  photometric  measurements  presented 
here  are  likely  contaminated  by  the  surrounding  star-forming 
regions  west  of  the  ULX  (see  the  HST  image  in  Eigure  3 
and  the  OM  image  of  Eigure  4).  As  a  result,  in  constructing 
the  SED  we  assume  these  are  upper  bounds  to  the  ULX 
emission. 

X-ray  to  UV  band.  The  UV  band  is  responsible  for  producing 
the  IR  emission  lines,  however,  it  is  also  unfortunately  the  most 
ambiguous  to  fit.  We  found  that  both  X-ray  models  that  have 
power-law  components  extending  to  lower  energies  (PLMCD 
and  BPL),  overpredict  the  measured  fiuxes  if  extrapolated  into 
UV  and  optical.  Therefore,  we  use  the  Modified  PLMCD  model 
to  extrapolate  the  X-ray  spectrum  to  lower  energies.  Here,  the 
disk  emission  can  be  approximated  as  a  simple  power  law  with 
photon  index  2/3  (or  energy  index  —1/3;  see  Mitsuda,  et  al. 
1984),  since  this  extrapolated  disk  component  is  a  superposition 
of  numerous  black  bodies  representing  the  integrated  emission 
from  the  disk  at  the  various  temperatures  and  radii.  At  optical 
energies,  this  model  predicts  a  drop  in  flux  due  to  the  finite  size 
of  the  accretion  disk  (see  Eigure  5). 

Optical  band.  The  disk  model  was  further  extrapolated  into 
the  optical,  where  the  companion  star  dominates  the  SED. 
KWZ  found  that  stars  within  the  spectral  type  range  B3Ib  to 
04V  were  consistent  with  optical  data.  Based  also  on  the  UV 
fluxes  presented  here,  we  chose  a  B2Ib  supergiant  as  the  optical 
counterpart.  The  UV  fluxes  seem  to  be  too  low  for  an  05V 
star.  The  model  for  the  B2Ib  star  is  part  of  the  ATLAS 9  Model 
Atmospheres  (Castelli  &  Kurucz  2004),  with  a  temperature  of 
18,500  K  and  a  luminosity  of  5  x  10^^  erg  s“^  Eigure  6  shows 
that  the  supergiant  emits  very  few  photons  that  can  impact  the 
ionization.  On  the  other  hand,  an  05  V  star  provides  a  significant 
flux  above  13.6  eV,  but  is  equally  deficient  above  the  He  ii  edge. 

Infrared  and  radio  band:  At  longer  wavelengths,  from  IR 
to  radio,  we  used  a  BPL  to  approximate  the  SED.  These 
photometric  data  are  viewed  as  upper  bounds  since  the  ULX  flux 
is  likely  contaminated  by  the  surrounding  star-forming  region. 
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Figure  6.  Representation  of  SED  constructed  in  Section  4.  The  radio  data  are  from  Miller  et  al.  (2005),  the  V-band  magnitude  for  the  optical  counterpart  of  the  ULX 
is  quoted  from  KWZ.  Our  measured  fluxes  from  IRAC,  MIPS,  OM,  and  GALEX  are  upper  bounds,  and  therefore  our  model  is  below  these  data  points.  The  two 
stellar  spectra  correspond  to  the  spectral  types  range  consistent  with  the  colors  and  magnitudes  in  KWZ.  The  X-ray  data  are  from  the  long  (~100  ks)  XMM-Newton 
observation  ID  0200470I0I.  The  Modifled  PLMCD  model  was  obtained  by  fitting  the  MCD  and  power-law  components  separately  as  explained  in  the  text.  The  other 
two  X-ray  models,  PLMCD  and  BPL,  are  also  shown  extrapolated  to  lower  energies,  they  clearly  overpredict  the  UV  and  optical  fluxes.  The  [O  iv]  edge  at  54.93  eV 
is  shown  as  a  vertical  line. 


5.  DISCUSSION 

5.7.  The  SED  of  the  ULX  in  Holmberg  II 

In  Section  4  and  Figure  6,  we  constructed  the  full  SED  for  the 
Holmberg  II  ULX  using  multi- wavelength  data  and  previously 
published  data.  The  bolometric  luminosity  for  this  SED  is  1 .34  x 
1040  g-i^  when  the  far-UV  spectrum  is  extrapolated  from 
the  X-rays,  using  the  Modified  PLMCD  model  fit.  Table  2 
illustrates  the  degeneracy  resulting  from  the  three  methods  for 
X-ray  spectral  fitting.  All  three  show  a  fairly  good  fit  to  the 
X-ray  data,  making  it  very  difficult  to  distinguish  between  the 
different  physical  interpretations  that  each  model  represents. 

The  optical  and  UV  fiuxes  plotted  in  Figure  6  can,  in  principle, 
come  from  three  different  sources:  the  stellar  companion,  the 
accretion  disk,  and  from  a  jet.  However,  the  overall  SED  shape 
and  in  particular  the  fact  that  the  optical  fiux  is  much  lower 
than  the  UV  fiuxes,  disfavor  the  jet  scenario  as  discussed  below. 
On  the  other  hand,  the  emission  from  the  disk  in  the  Modified 
PLMCD  model  is  too  faint  at  optical  and  UV  wavelengths  if 
extrapolated  from  X-rays  and  requires  an  additional  source. 
Previously  published  optical  data  and  our  UV  measurements 
suggest  that  a  ~B2Ib  star  is  the  source  of  these  fiuxes  as  the 
companion  to  the  accreting  BH.  We  explore  this  further  in 
Paper  II. 

The  fact  that  the  PLMCD  model  overpredicts  the  UV  and 
optical  fiuxes  if  extrapolated  to  these  wavelengths  is  expected, 
because  the  power  law  cannot  extend  too  much  below  the  peak 
of  the  accretion  disk  (e.g..  Done  &  Kubota  2006).  On  the  other 
hand,  it  is  intriguing  that  the  simple  BPL  model  provides  such 
a  good  fit  to  the  high-quality  X-ray  data.  Moreover,  the  location 
of  the  break  (~1  keV)  in  this  case  is  in  agreement  with  the 
jet  model.  Kaaret  et  al.  (2003)  fit  a  very  similar  model  to 
another  famous  ULX  in  NGC  5408.  It  has  been  shown  that 


jets  can  contribute  significantly  to  the  X-ray  emission,  both 
by  optically  thin  synchrotron  emission  and  Compton  scattering 
of  the  photons  from  the  disk  (or  the  companion)  by  relativistic 
electrons  in  the  jet  (Georganopoulos  et  al.  2002).  However,  even 
if  the  power  law  is  less  steep  compared  to  the  PLMCD  model,  it 
is  still  overpredicting  the  UV  and  optical  fiuxes  and  is  therefore 
not  the  preferred  scenario. 

Much  better  constraints  are  provided  by  photoionization 
modeling  of  the  Spitzer  spectra.  The  [O  iv]  line  is  only  produced 
by  high-energy  photons  (^54  eV),  and  is  clearly  associated 
with  the  ULX  (Figure  3).  For  instance,  the  “classical”  PLMCD 
model  predicts  a  power-law  component  with  photon  index  of 
2.4  dominating  into  the  UV.  This  model  has  the  strongest 
emission  of  the  three  models  in  the  far-UV  and  soft  X-rays, 
and  should  produce  the  strongest  [O  iv]  emission.  In  this  case, 
photoionization  modeling  can  be  used  to  predict  the  [O  iv]  fiux 
predicted  from  such  a  model,  and  this  predicted  fiux  can  be 
compared  with  our  observed  values.  In  a  companion  paper 
(Paper  II),  we  present  the  results  of  such  a  photoionization 
analysis  in  an  attempt  to  distinguish  between  the  three  X-ray 
models  presented  here. 

5.2.  Characterization  of  the  Optical/Mid-IR  Line  Ratios 

The  ionization  potential  of  [Oiv]  is  just  above  the  Hen 
edge  so  both  likely  result  from  similar  regions  in  the  nebula 
and  a  quantitative  comparison  between  previously  published 
optical  observations  and  the  mid-IR  observations  can  be  made. 
KWZ  found  that  the  He  ii  recombination  line  is  consistent  with 
photoionization  by  the  ULX.  The  fiux  of  our  [O  iv]  line  is  almost 
three  times  higher  than  the  He  ii  fiux  measured  by  KWZ  (2.7  x 

1036  gj.g 

We  compared  the  mid-IR  line  ratios  of  the  ULX  with  those 
in  AGN.  NGC  4395  is  one  of  the  few  known  galaxies  that 
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Figure  7.  Ratio  line  diagnostics  comparison  between  Holmberg  II ULX  and  other  extragalactic  ionizing  sources:  AGN,  starburts,  etc.  Its  position  in  both  plots  indicates 
a  source  very  similar  to  AGNs  and  LINERs.  (a)  Diagnostic  that  uses  the  [O  iv]  line  normalized  to  the  [Ne  ii]  line.  The  AGN,  starburts,  and  ULIRG  data  are  taken  from 
Genzel  et  al.  (1998).  Many  of  the  starburts  and  ULIRGs  are  actually  upper  limits  (see  Figure  3  of  Genzel  et  al.  1998).  For  Holmberg  II UFX,  we  show  the  lower  limit 
because  [Neii]  was  not  detected,  (b)  Diagnostic  plot  that  uses  the  [S  iii]  and  [Sin]  lines  to  distinguish  between  different  sources.  We  plot  all  the  data  presented  by 
Dale  et  al.  (2006,  their  Figure  5),  which  include:  Seyferts,  LINERs,  starburt  galaxies,  Hii  regions,  and  extranuclear  star-forming  regions.  For  Holmberg  II  ULX,  we 
show  the  lower  and  upper  limits  because  [Neii]  and  [S  in],  respectively,  were  not  detected. 


is  likely  powered  by  an  IMBH  (Eilippenko  &  Ho  2003).  It 
has  a  bolometric  luminosity  very  close  to  our  ULX,  and  high- 
excitation  emission  lines  were  detected.  The  narrow-line  He  ii 
measured  by  Kraemer  et  al.  (1999)  for  NGC  4395  relative  to 
H)0  is  0.12,  a  value  typical  for  Seyfertl  galaxies  (Dietrich  et  al. 
2005).  Eor  the  Holmberg  II  ULX,  Heii/H^  ^  0.2  (KWZ), 
slightly  higher  than  for  NGC  4395.  The  [Oiv]  measured  from 
NGC  4395  is  almost  16  times  larger  than  the  He  ii  flux  (Kraemer 
et  al.  1999).  However,  the  gas  density  is  close  to  the  critical 
density  for  this  line  (10"^  cm“^) — much  higher  than  the  density 
of  the  gas  surrounding  the  Holmberg  II  ULX  (KWZ  estimate 
10  cm“^,  based  on  the  line  surface  brightness  predicted  with 
CLOUDY).  The  [Oiv]  is  expected  to  increase  as  it  approaches 
the  critical  density,  while  the  Heii  is  not  sensitive  to  the  gas 
density.  An  even  more  important  difference  is  the  abundance 
of  O,  which  is  ~5  times  higher  for  NGC  4395  (Kraemer 
et  al.  1999).  The  dependence  of  the  [Oiv]  on  the  abundance  is 
expected  to  be  roughly  linear.  In  spite  of  these  two  differences, 
we  consider  the  [O  iv]  (and  He  ii)  detected  from  the  Holmberg  II 
ULX  to  be  consistent  with  photoionization  by  a  high  energy 
source  similar  to  the  AGN  in  NGC  4395. 

The  infrared  [Oiv]  line  (together  with  [Ne  v])  is  a  signature 
of  high  excitation  usually  associated  with  actively  accreting 
black  holes  (e.g..  Lutz  et  al.  1998).  This  line  is  often  used  to 
disentangle  the  starburst  emission  in  composite  type  galaxies, 
and  is  a  good  indicator  of  the  AGN  power  (Genzel  et  al.  1998; 
Sturm  et  al.  2002;  Satyapal  et  al.  2004,  2007;  Smith  et  al. 
2004;  Melendez  et  al.  2008).  Genzel  et  al.  (1998)  and  Sturm 
et  al.  (2002)  have  used  the  ratio  of  [Oiv]/[Neii]  to  distinguish 
between  starbursts  and  AGN-dominated  galaxies  in  the  mid- 
IR.  In  Ligure  7(a),  we  compare  the  Holmberg  II  ULX  with 
data  published  by  Genzel  et  al.  (1998).  The  comparison  clearly 
indicates  high  excitation  for  the  ULX,  similar  to  AGN. 

Other  infrared  line  diagnostics  have  been  proposed  by  Dale 
et  al.  (2006,  2009)  to  disentangle  AGN  type  ionizing  sources 
from  star  formation,  using  the  ratios  [Siii]  33.48  /xm/[Siii] 


34.82  /xm  and  [Neiii]/[Neii].  Ligure  7(b)  shows  the  four 
regions  deflned  by  Dale  et  al.  (2006)  in  their  Ligure  5,  and  their 
published  data  for  different  sources  in  the  SINGS  sample.  The 
Holmberg  II  ULX  is  located  well  within  region  I,  together  with 
most  AGNs  and  LINERs.  We  note  that  Vazquez  et  al.  (2007) 
found  similar  results  in  their  Spitzer  IRS  study  of  the  NGC 
4485/90  ULXs.  However,  they  did  not  And  clear  evidence  for 
[O  iv]  in  their  sample  of  objects.  In  examining  this  data  further, 
we  And  slightly  discrepant  ratios  from  Vazquez  et  al.  (2007). 
Lrom  our  analysis,  only  ULX-2  in  their  sample  actually  falls 
marginally  inside  region  I,  the  other  flve  ULXs  being  located 
in  regions  II  and  III.  In  communicating  with  the  authors,  we 
And  that  the  discrepancy  is  likely  do  to  differences  in  reduction 
pipeline  used  for  processing  or  the  apertures  used  in  either 
analysis. 

Both  diagnostic  diagrams  suggest  that  the  [O  iv]  emission  in 
Holmberg  II  is  dominated  by  the  ULX  and  that  the  surrounding 
star-forming  regions  provide  little  contamination  to  this  line. 
While  this  comparison  provides  good  evidence  that  the  line 
ratios  for  the  ULX  are  very  similar  to  AGN,  we  note  that  AGNs 
are  different  than  ULXs,  both  in  their  intrinsic  emission  and 
their  environment.  Lor  example,  AGNs  emit  copiously  in  the 
UV  and  optical,  and  are  surrounded  by  gas  at  high  densities. 

Dasyra  et  al.  (2008)  found  a  good  correlation  between  the 
luminosity  of  the  [Oiv]  line  and  the  mass  of  the  central  BH 
in  standard  AGN.  We  estimate  a  BH  mass  of  3.5  x  Mq 
if  this  correlation  is  extrapolated  to  the  Holmberg  II  ULX.  We 
can  also  estimate  the  BH  mass  using  the  inner  disk  temperature 
from  the  Modifled  PLMCD  model  (0.38  keV).  We  obtain  a  BH 
mass  of  994  Mq  .  However,  the  flux  of  the  MCD  disk  component 
does  not  dominate  the  X-ray  emission  (is  only  49%  of  the  total 
flux).  This  scenario  would  imply  that  the  ULX  is  not  in  a  typical 
high  state,  and  therefore  the  inner  disk  temperature  cannot  be 
used  reliably  to  estimate  the  mass  of  the  BH  (Stobbart  et  al. 
2006).  Linally,  using  the  bolometric  luminosity  calculated  from 
our  SED  (1.34  x  10^^  erg  s“^),  we  estimate  that  the  lower  limit 
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to  the  mass  of  the  central  BH,  within  the  Eddington  regime,  is 
85  Mq. 

6.  CONCLUSIONS 

The  detection  of  the  [O  iv]  line  from  the  Holmberg  II  ULX 
indicates  strong  UV  and  soft  X-ray  emission.  This  line  is 
commonly  used  to  quantify  the  AGN  power  in  Seyfert  galaxies. 
The  ratios  of  this  line  to  lower-ionization  lines  are  comparable 
to  those  in  AGNs  and  suggest  that  the  [Oiv]  emission  in 
Holmberg  II  results  from  the  ULX  and  not  from  the  surrounding 
star- forming  regions.  The  mid-IR  line  ratios  indicate  that  the 
intrinsic  luminosity  of  the  ULX  must  be  very  high. 

The  SED  constructed  based  on  the  X-ray  spectral  modeling, 
UV  and  infrared  photometry,  and  previously  published  radio  and 
optical  data  suggest  that  the  ULX  companion  is  likely  similar  to  a 
B2Ib  supergiant.  The  accretion  disk  is  not  expected  to  contribute 
very  much  in  the  optical.  The  overall  shape  of  the  SED  seems  to 
favor  the  accretion  disk  emission  as  the  dominant  component  in 
the  soft  X-ray  and  far-UV,  and  is  therefore  likely  responsible  for 
the  photoionization.  In  the  accompanying  paper  (Paper  II),  we 
perform  detailed  photoionization  modeling  to  provide  further 
evidence  that  the  detected  [Oiv]  line  is  consistent  with  the 
ionizing  emission  from  the  ULX,  and  to  better  constrain  the 
three  X-ray  models  presented  here. 
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